Supersolid hardcore bosons on the triangular lattice 
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We determine the phase diagram of hardcore bosons on a triangular lattice with nearest neighbor 
repulsion, paying special attention to the stability of the supersolid phase. Similar to the same model 
on a square lattice we find that for densities p<l/3orp>2/3a supersolid phase is unstable and 
the transition between a commensurate solid and the superfluid is of first order. At intermediate 
fillings l/3<p<2/3we find an extended supersolid phase even at half filling p = 1/2. 
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Next to the widely observed superfluid and Bose- 
condensed phases with broken U{\) symmetry and "crys- 
talline" density wave ordered phases with broken trans- 
lational symmetry, the supersolid phase, breaking both 
the U{\) symmetry and translational symmetry has been 
a widely discussed phase that is hard to find both in 
experiments and in theoretical models. Experimentally, 
evidence for a possible supersolid phase in bulk ^He has 
recently been presented |l| , but the question of whether 
a true supersolid has been observed is far from being set- 
tled 0, 01 J leaving the old question of supersolid behavior 
in translation invariant systems 0, 0] unsettled for now. 

More precise statements for a supersolid phase can be 
made for bosons on regular lattices. It has been pro- 
posed that such bosonic lattice models can be realized 
by loading ultracold bosonic atoms into an optical lattice, 
where the required longer range interaction between the 
bosons could be induced by using the dipolar interaction 
in chromium condensates |a| , or an interaction mediated 
by fermionic atoms in a mixture of bosonic and fermionic 
atoms . With the recent realization of a Bose-Einstein 
condensate (BEG) in Chromium atoms these exper- 
iments have now become feasible, raising the interest in 
phase diagrams of lattice boson model, and particularly 
in the stability of supersolids on lattices. 

The question if a supersolid phase is a stable ther- 
modynamic phase for lattice boson models has been 
controversial for many years. Analytical calculations 
using mean-field and renormalization group methods 
have predicted supersolid phases for many 
models, including for the simplest model of hardcore 
bosons with nearest neighbor repulsion on a square lat- 
tice with Hamiltonian 



H 



(1) 



where a| (a^) creates (destroys) a particle on site z, 
t denotes the nearest-neighbor hopping, V a nearest- 
neighbor repulsion, and fi the chemical potential. Sub- 
sequent numerical investigations using exact diagonal- 
ization and quantum Monte Carlo (QMC) algorithms 
[II ElIIl EE 113 have shown that for this model, the 
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FIG. 1: Zero-temperature phase diagram of hardcore bosons 
on the triangular lattice in the canonical ensemble obtained 
from quantum Monte Carlo simulations. The regions of phase 
separation are denoted by PS. The insets exhibit the density 
distribution inside the solid phases for p = 1/3 (lower panel), 
and p = 2/3 (upper panel). 



supersolid phase is unstable and phase separates into su- 
perfluid and solid domains at a first order (quantum) 
phase transition. Recently, this occurrence of a first or- 
der phase transition was explained by showing that a 
uniform supersolid phase in a hardcore boson model is 
unstable towards the introduction of domain walls, low- 
ering the kinetic energy of the system by enhancing the 
mobility of the bosons on the domain wall In a re- 
lated work it has been proposed that superfluid domain 
walls might be an explanation for the experimental ob- 
servation of possible supersolidity in Helium p, 18] . 

To stabilize a supersolid on the square lattice, the ki- 
netic energy of the bosons in the supersolid has to be 
enhanced either by sufficientl y r educing the on-site in- 
teraction to be less than 4:V [l3|, by adding additional 
next-nearest-neighbor hopping terms [l^ , or by forming 
striped solid phases with additional longer-ranged repul- 
sions [ll[l3- 

In this Letter we will consider the interplay of super- 
solidity and frustration by studying the hardcore boson 
model on a triangular lattice. In the classical limit 
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FIG. 2: Zero-temperature phase diagram of hardcore bosons 
on the triangular lattice in the grand canonical ensemble ob- 
tained from quantum Monte Carlo simulations. Second order 
phase transitions are denoted by solid lines, whereas first- 
order transitions are denoted by dashed lines. The system is 
half-filled for p^jV = 3. 
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FIG. 3: Density of hardcore bosons on the triangular lattice 
as a function of {i along lines of constant values of tjy . The 
inset displays the jump in the density as a function of t for 
[ijV = 4 at tjV ~ 0.165. Only densities p > 1/2 are shown 
since the phase diagram is symmetric around half filling. 



t — two solid phases exist at fillings p = 1/3 (and 
p — 2/3), where one of three sites is filled (empty) in a 
\/3 X \/3 ordering with wave vector Q = (47r/3, 0) [20j . 
shown in the insets of Fig. ^ At half filling {p = 1/2), 
where the square lattice shows a solid ordering with wave 
vector (tt, tt), the solid order is frustrated on the trian- 
gular lattice, and the classical model has a hugely degen- 
erate ground state with an extensive zero-temperature 
entropy [2]| . 

The question arises whether this degeneracy of the 
classical system at half filling is lifted when quantum dy- 
namics is added at a finite hopping parameter t, and 
which phase gets stabilized. Mean-field studies have pre- 
dicted a supersolid phase Given the questionable 
reliability of mean-field calculations in the case of the 
square lattice model a numerical check is needed. Indeed 
Green's function Monte Carlo (GFMC) simulations on 
small lattices 23] have indicated the absence of a super- 
solid phase at half filling, but again on such small lattice 
numerical results can also be misleading as in the square 
lattice model [l3| . 

We have thus performed a series of high-accuracy nu- 
merical QMC calculations on large lattices using stochas- 
tic series expansions with global directed-loop up- 
dates [i^ for the hardcore boson model on the triangular 
lattice and show the phase diagram in Fig. ^ and Fig. [3 
for the canonical and grand-canonical ensemble, respec- 
tively. The main results are that for fillings p < 1/3 and 
p > 2/3 a supersolid is unstable towards phase separation 
by exactly the same domain-wall proliferation mechanism 
through which the square lattice supersolid is unstable at 
all fillings p ^ 1/2. In contrast, for intermediate densi- 
ties l/3<p<2/3we find that the degeneracy of the 
frustrated classical model is indeed lifted and a stable 



supersolid phase emerges. The phase diagram in Fig. [21 
is similar to the mean- field phase diagram [i^l, albeit 
with a substantially reduced supersolid region. The su- 
persolid is stable even at half filling, contradicting the 
small-lattice GFMC results of Ref. 23. 

We will now discuss the phase diagrams in more de- 
tail, starting with simple limits. Considering the sin- 
gle boson (hole) problem, one can show that the lattice 
is empty for p < po — — 6i and completely filled for 
p > Pi = 6{t + V). For large values of t/V, the bosons 
are superfluid, with a finite value of the superfluid density 
PS, which we measure through the winding number fluc- 
tuations W of the world fines ^ a.s ps = (W^) / {4/31) . 
Two solid phases emerge upon lowering t/V with rational 
fiUings 1/3 and 2/3, respectively. Both are characterized 
by a finite value of the density structure factor per site, 
Siq)/N = (pqpjj), where p^ = {l/N)J2,niexp{iqr,) 
at wave vectors ±Q = ±(47r/3, 0), corresponding to 
the VS X \/3 ordering wave vector. The maximum ex- 
tent of the solid phases is reduced by quantum fluctua- 
tions from the mean-field value of {t/V)c — 0.5 down to 
(t/F)c ==0.195 ± 0.025. 

Since the phase diagram is symmetric when inter- 
changing particles with holes (p — > 1 — p) we restrict 
our discussion from now on to p > 1/2 and plot the den- 
sity p as a function of chemical potential p for cuts at 
constant t/V in Fig.|Hl For t/V — 0.1 we clearly observe 
a plateaux corresponding to the p = 2/3 (p = 1/3) phase 
with broken translational symmetry. The approach to 
this plateaux from p<2/3(p>l/3)is continuous, indi- 
cating a second order phase transition, while for p > 2/3 
(p < 1/3) we see a jump caused by a first order phase 
transition. Measuring the density structure factor S'(q) 
and the superfluid density in Fig. 01 we identify this as a 
first order phase transition between the solid and super- 
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FIG. 4: Static structure factor S'(Q) for hardcore bosons 
on the triangular lattice as a function of fi along a line of 
constant t/V = 0.1. The inset shows the behavior of the 
superfluid density ps- 





FIG. 5: The p — 2/3 solid doped with bosons, a) additional 
bosons (open circles) added on top of the solid, b) lining the 
bosons up costs no additional potential energy, c) shifting the 
lower half of the lattice introduces a domain wall (dashed line) 
at no cost, but now d) the additional particles can hop freely 
across the domain wall, gaining additional kinetic energy. 
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FIG. 6; Finite size scaling behavior of the static structure 
factor ^(Q) and the superfluid density ps for hardcore bosons 
on the triangular lattice at t/V — 0.1 and half filling (p = 1/2, 
p/V — 3). Dashed lines indicate extrapolations to the infinite 
lattice. 
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FIG. 7: Static structure factor S(Q) for hardcore bosons on 
the triangular lattice as a function of t at half filling {p — 1/2 
and p/V = 3.) The inset shows the behavior of the superfluid 
density ps and the kink at t/V ^ 0.12, indicated by an arrow. 



fluid phases. 

The situation here is the same as in the square lattice 
model, where doping the solid leads to phase separation 
at a first order phase transition. The strict arguments 
for instability of a supersolid phase in the square lattice 
can also be applied here: the uniform supersolid is 
unstable towards the introduction of domain walls as we 
illustrate in Fig. El We start by adding L/3 additional 
bosons to the solid at density p = 2/3 [Fig. |3Ji)], which 
corresponds to an infinitesimal density in the thermody- 
namic limit. These bosons can gain a kinetic energy of 
—6t^/V per boson by second order hopping processes. 
Placing these additional bosons along a line, as shown in 
Fig. [Sja) costs no additional potential energy, and we can 
even shift one half of the lattice by one lattice spacing, in- 



troducing a domain wall as shown in Fig. (SJ:), again at no 
cost in potential energy. But now, the additional bosons 
can gain kinetic energy of —t per boson by hopping freely 
across the domain wall, which lowers the energy of the 
domain wall state compared to the bulk supersolid, and 
hence the supersolid phase is unstable. 

A different situation exists for p < 2/3, since there is 
no symmetry around p — 2/3. Here, forming a domain 
wall would cost extra potential energy, and a supersolid 
phase can thus be stabilized. To demonstrate the exis- 
tence of this supersolid even at half filling, we show the 
finite size scaling of ps and <S'(Q) in Fig. (HJ both of which 
extrapolate to finite values. Intervening the solid phases 
atl/3<p<2/3we hence find an extended supersolid 
phase, where both the superfluid density and the density 
structure factor take on flnite values. Fig. \7\ shows ps 
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in ultra-cold atoms on optical lattices. 
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FIG. 8: Static structure factor S'(Q) for hardcore bosons on 
the triangular lattice as a function of t along lines of constant 
^/V — 3.4 and fi/V = 4. The inset shows the superfluid 
density ps, exhibiting a kink at t/V ~ 0.125 for /i/V = 3.4. 



and S'(Q) as functions of t/V at half filling, indicating 
a continuous quantum phase transition from the super- 
solid to the superfluid at t/V ~ 0.115. We observe a 
kink in psit) near the transition point, marked by an ar- 
row in Fig. Away from half-filling, the extend of the 
supersolid phase slightly increases, as shown in Fig. |21 
Moreover, the kink in ps{t) at the supersolid-superfluid 
transition becomes more pronounced, being clearly visi- 
ble for p/V = 3.4 in Fig.|Hl Eventually, for p/V > 3.95, 
the supersolid phase ceases to be stable, giving rise to 
a direct first order transition between the solid and the 
superfluid. This is reflected in discontinuities of both ps 
and S{Q) in Fig. [SJ as well as in the density p [Fig. 

To summarize, we have demonstrated that, in contrast 
to square lattice hardcore boson models, an extended su- 
persolid phase exists on the triangular lattice without the 
need for longer-range or softcore interactions, albeit in a 
smaller region than predicted by mean-field calculations 
|22l | and partially contradicting previous simulations on 
smaller lattices . This supersolid phase in the density 
regime 1/2 < p < 2/3 emerges from the hugely degen- 
erate disordered ground state of the frustrated classical 
model (in the t — limit) when the quantum mechani- 
cal hopping is turned on. This illustrates an intriguing 
mechanism by which a quantum system can avoid frus- 
tration: while A^/3 of the bosons, on an A^-site lattice 
form a non- frustrated solid at wave vector (47r/3,0) and 
break translational symmetry, the remaining N{p~ 1/3) 
bosons delocalize and break the U{1) gauge symmetry, 
forming a superfluid Bose-condensate on top of the solid 
with density p — 1/3, thus realizing a supersolid phase. 

Since, in contrast to the square lattice, the triangular 
lattice model does not need additional longer-ranged re- 
pulsion or hopping terms, nor a reduction of the on-site 
interaction |l7j , the triangular lattice might be preferred 
over the square lattice when looking for supersolid phases 
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